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a  turbulent  velocity  probe,  and  a  wave-height  probe  in  order  to  investigate 
the  statistics  and  environmental  causes  of  amplitude  modulation.  The  temporal 
variations  of  the  sound  amplitude  were  studied  during  four  twenty-minute  runs 
over  a  range  of  2  meters  at  varying  depths  from  near  the  surface  to  l4  meters 
during  sea-state  one  conditions  in  water  of  16  meters-depth.  Analysis  in¬ 
dicated  mean  percentage  amplitude  modulation  of  approximately  5%  with  a 
variance  of  approximately  10”3.  Temporal  correlation  times  were  about 
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INTRODUCTION 


A.  ACOUSTIC  AMPLITUDE  MODULATION  IN  A  RANDOM  MEDIUM 

Consider  an  underwater  acoustic  wave  propagated  at  a  constant 
frequency  and  a  fixed  amplitude  in  isotropically  homogeneous,  un¬ 
disturbed  water  with  no  scatterers.  In  the  far  field  of  the  trans¬ 
ducer,  the  wave  will  arrive  at  the  receiver  unchanged  except  for  a 
diminution  due  to  spherical  spreading  loss  and  absorbtion.  The 
amplitude  will  not  vary  about  the  peak-to-pcak  value  received. 

For  low  values  of  absorbtion,  this  type  of  acoustic  wave  is  ade¬ 
quately  described  by  the  wave  equation: 

V  ZP  ~  _L 

I  Cx?l*- 

Near  the  ocean  surface,  a  wave  will  encounter  random  inhomo- 
genieties  which  cause  a  varying  speed  of  sound  resulting  in  re¬ 
fraction.  The  wave  equation  must  therefore  be  modified  to  include 
temporal  and  spatial  changes  in  the  speed  of  sound.  Such  changes 
may  be  exoiessed  as  a  change  in  the  index  of  refraction  U'r,t), 
where  U(r,t)  is  the  ratio  of  average  sound  speed  to  local  sound 
speed:  C 


ju(r.i)  r  — 1— 


which  may  be  expressed  as: 

ytff.t )  =  I  *  ocntr.t) 
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where  Of  is  the  root-mean-square  variation  of  the  index  of  re- 
fraction  from  unity,  and  n(r,t)  is  the  normalized  variation  of  [JL 


from  unity.  The  quantity  ft  has  been  shown  to  be  a  small  number, 
equal  to  1.18  x  10  ^  in  laboratory  experiments  using  heated  water 
(Stone  and  Mintzer ,  1965)  and  to  2.55  x  10  ^  in  ocean  experiments. 
(Whitmarsh,  Skudrzyk  and  Urick,  1957.) 

The  effect  of  the  random  variation  of  the  index  of  refraction 
on  an  acoustic  wave  may  be  shown  by  a  variation  in  the  amplitude 
of  a  received  signal,  noticable  as  a  low  frequency  modulation  of 
the  carrier  frequency.  Such  modulation  may  be  expressed  as  a 
fraction  of  the  maximum  undisturbed  amplitude  of  the  received 
signal  as  shown  below: 

instantaneous  fractional  modulation  = 

wher  e 

VMAX  =  the  unmodulated  voltage  value 

VMOD  =  the  instantaneous  modulated  voltage  value 

The  inhomogeniet  ies  which  cause  the  variations  in  the  re¬ 
fractive  index,  and  ultimately  result  in  amplitude  modulation,  are 
temporal  and  spatial  variations  in  temperature  and  salinity  micro- 
structure,  random  bubble  populations,  particulate  matter,  varying 
forms  of  oceanic  life  and  those  inhomogen iet ies  caused  by  man. 

Each  factor  will  vary  with  location  and  ambient  conditions  in  the 
ocean.  The  effect  on  sound  transmission  may  also  be  a  function  of 
sound  frequency. 

Particle  velocit  in  the  ocean  has  a  great  effect  on  the  above 
parameters,  generally  in  determining  their  temporal  and  spatial 


VMAX  -  VMOD 
VMAX 
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distribution.  There  are  two  main  types  of  oceanic  particle  velocity, 
that  due  to  convective  turbulence  and  that  due  to  the  effect  of  sur¬ 
face  waves.  Convective  turbulence  is  fairly  slow  and  affects  the 
distribution  of  temperature  and  salinity.  Such  concepts  as  tem¬ 
perature  "patches"  and  lenticular  volumes  of  near  unitemperature 
are  due  to  convective  turbulence. 

Turbulence  arising  from  the  effect  of  surface  wave  action  is 
generally  responsible  for  random  variations  in  populations  of 
scatterers  and  for  small  scale  variations  in  temperature  and 
salinity  distribution. 

Distr ibtuions  and  variations  of  all  parameters  affected  by 
wave  turbulence  would  show  a  lessening  effect  with  increasing 
depth  due  to  the  e  '  decay  of  wave  amplitude.  Convective  tur¬ 
bulence  could  exist  at  all  depths  near  the  surface  and  its  de¬ 
pendence  on  depth  is  unknown.  These  mechanisms  are  always  present 
near  the  surface  of  the  ocean  and  are  a  principal  factor  in  the 
modulation  of  sound.  They  also  lead  to  problems  in  analysis.  For 
example,  a  stationary  temperature  patch  might  show  temporal  char¬ 
acteristics  due  to  displacement  across  a  sensor  with  wave  action. 

Statistical  analysis  of  the  amplitude  modulation  occurring  in 
inhomogeneous  transmission  has  centered  about  the  coefficient  of 
variation,  CV,  defined  by  Stone  and  Mintzer  (1961): 

(cv/.  2 

tc ;  <p>1 

where  CV  =  the  coefficient  of  variation 

P  -  the  amplitude  of  the  received  signal 
O  =  the  ensemble  average. 
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The  same  coefficient  of  variation  is  defined  equivalently  by 


i 


I 


| 

| 

I 


Jl 

■ 


I; 


I 


Campanella  and  Favret  (1969)  as: 


(cvy 


SSU)  -  <&u)» 


<SU)> 
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where  S(t)  is  the  amplitude  of  the  received  signal. 

The  relationship  between  the  coefficient  of  variation  as  de¬ 
fined  by  Stone  and  Mintzer ,  (CV)  ,  and  the  fractional  modulation  as 

defined  in  this  thesis  is  shown  in  Appendix  I. 

2 

The  relation  between  (CV)  and  the  RMS  variation  of  the  index 
of  refraction  O',  the  range  r  and  the  size  of  the  inhomogen ieties 
has  been  variously  described,  as  by  Stone  and  Mintzer  (1965)  de¬ 
pending  on  the  (presumed)  shape  of  the  spatial  correlation  function 
of  the  index  of  refraction  R (p)  where: 

=  <n(r)  n(r*-p)> 

If  R(p)  is  Gaussian, 

*R(p>  =  e*  p  (  -  P  / a.a  ) 


then  the  coefficient  of  variation  may  be  derived  for  the  conditions 

kr  »  1 
ka  »  1 


that 


and,  r/ka  «  1 


as : 


r 


where  k 

o 

a. 

a 

r 


space  wave  number  of  the  sound  frequency 
RMS  variation  of  the  index  of  refraction 
radius  of  the  scattering  inhomogeniety 
transducer  to  receiver  range. 
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If  R(p)  is  exponential; 


R  (?)  =  exp  (- /°/ol,) 

2 

then  the  numerical  factor  in  the  expression  for  (CV)  above  would 
be  unity.  Although  an  exponential  correlation  function  would  not 
satisfy  the  requirements  of  non-zero  slope  at  p  =  0,  the  experi¬ 
mental  work  (Stone  and  Mintzer ,  1961  and  1965,  and  Campanella  and 
Favret,  1969)  seems  to  indicate  an  exponential  correlation 
function.  The  dilemna  cai.  be  resolved  by  using  a  modified  ex¬ 
ponential  function  to  satisfy  the  p  =  0  condition. 

The  temporal  correlation  function  for  the  variation  of  the 
index  of  refraction, 

M(*r}  =  >0  Ct t  T) y 

has  been  compared  to  the  temporal  correlation  function  for  the 
variation  of  sound  pressure  amplitudes, 

and  the  relation: 

M  (t )  -  G>(t) 

which  was  derived  theoretically  by  Mintzer  (1954)  and  was  shown 
to  exist  experimentally  by  Campanella  and  Favret  (1969)  in  a  con¬ 
trolled  laboratory  environment  using  convective  temperature 
differences  to  generate  a  varying  index  of  refraction. 

Experimental  work  in  the  ocean  environment  has  included  pro¬ 
pagation  at  frequencies  from  800  Hz  (R.M.  Kennedy,  1969)  and 
60  kHz  (Whitmarsh  et  al.,  1957)  at  distances  of  1000  yds 
(Whitmarsh  et  al.)  and  twenty-five  miles.  (Kennedy.) 
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Kennedy  found  the  probability  distribution  function  of  the 
logarithm  of  the  amplitude  fluctuations  to  be  normally  distributed, 
and  an  exponential  temporal  correlation  function  with  an  average 
decorrelation  time  to  1/e  in  12  minutes.  There  was  no  attempt  to 
take  measurements  of  the  environment  during  his  experiment. 
Whitmarsh,  Skudrzyk  and  Urick  found  that  the  fluctuation  in  amp¬ 
litude  tended  to  increase  with  range,  and  decried  a  lack  of  good 
information  concerning  the  microthermal  structure  of  the  ocean. 

B.  SCOPE  AND  INTENT  OF  THE  EXPERIMENT 

The  purpose  of  the  experiment  was  to  measure  and  analyze  the 
temporal  variation  in  the  fractional  modulation  of  an  acoustic  wave 
over  an  oceanographically  defined  volume  in  the  ocean.  The  tem¬ 
poral  characteristics  of  the  environmental  factors  of  temperature, 
salinity,  turbulent  velocity,  sound  speed  and  wave  height  were 
measured  simultaneously  with  fluctuations  in  sound  amplitude. 

A  simple  one-transducer  and  hydrophone  arrangement  over  a 
short  transmission  length  of  two  meters  was  adequate  to  measure 
variations  in  the  amplitude  of  the  transmitted  signal.  The  short 
range  was  desirable  because  of  the  need  for  system  rigidity  and  an 
adequate  definition  of  the  environmental  factors  occurring  during 
each  twenty-minute  run.  The  length  of  each  run  was  the  minimum 
time  necessary  to  statistically  define  the  low  frequency  environ¬ 
mental  factors  which  caused  the  amplitude  modulation. 

The  acoustical  and  oceanographic  analysis  at  each  depth  was 
concerned  with  the  statistical  and  temporal  characteristics.  The 
results  desired  were  the  mean  and  variance,  the  probability  density 
function,  the  temporal  correlation  and  the  power  spectral  density 


of  each  variable. 


11 


The  acoustic  results  will  later  be  cross-correlated  with  the 


varying  environmental  factors  to  determine  the  relative  effect  of 
each  parameter  upon  the  amplitude  modulation. 
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DESIGN  OF  THE  EXPERIMENT 


II  . 


A.  DETERMINATION  OF  THE  CONFIGURATION  OF  THE  EXPERIMENTAL  APPARATUS 


1 .  Conceptual  Design  of  the  Experiment 

In  order  to  measure  environmentally  caused  amplitude  modu¬ 
lation  over  a  clearly  defined  volume  in  the  ocean,  a  high  frequency 
acoustic  carrier  wave  was  directed  at  a  receiving  hydrophone 
mounted  upon  the  acoustic  axis  at  a  separation  of  two  meters.  The 
received  signal  was  processed  for  the  amplitude  modulation  occurring 
in  each  twenty-minute  run. 

In  support  of  the  acoustic  investigation,  an  electro¬ 
magnetic  flowmeter,  three  thermistors,  a  Bi ssett -Herman  apparatus 
with  a  salinometer  and  a  temperature  probe  and  a  Ramsey  probe  with 
a  velocimeter  and  a  wave-height  probe  were  mounted  about  the 
acoustic  apparatus  to  define  the  environment  during  each  run. 


2 .  The  Acoustic  Carrier  Wave 


The  carrier  wave  utilized  in  the  experiment  was  a  c 
tinuous-wave  60  kHz  constant  amplitude  signal.  This  frequenc 
selected  in  order  to  maximize  the  effect  of  random  bubble  popu¬ 
lations  which  were  anticipated  to  be  greatest  near  the  ocean  sur¬ 
face  . 

As  may  be  seen  from  Figure  1 1 - 1 ,  the  60  kHz  carrier  fre¬ 
quency  was  near  resonance  for  the  most  populous  bubble  sizes  near 
the  surface.  The  resonance  curves  were  calculated  on  the  basis  of 


the  formula: 
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70 U  Bubble  Radius 
60U  Bubble  Radius 


Depth  (meters) 


where 


=  resonance  frequency  { Hz ) 


f 

o 

R  =  radius  of  the  bubble  (cm) 
o 

Fq  =  ambient  pressure  at  the  bubble  depth,  equals  10^ 

2  62 

dynes/cm  at  the  surface  with  10  dynes/cm“  added 
for  each  10  meters  of  depth 

3 

=  density  of  water  (1.0  gram/cm'  ) 
y  -  ratio  of  specific  heats  (assumed  to  be  1.4) 
c{0l  *  =  correction  for  surface  tension  and  thermal 
conductivity  equals  I.O296  for  60  kHz 
3 .  Determination  of  the  Acoustic  Range 

The  most  stringent  requirements  for  the  distance  between 
transducer  and  hydrophone  arose  from  the  desirability  that 

(a)  the  hydrophone  be  in  the  iar  field  of  the  transducer. 

(b)  the  volume  of  water  under  investigation  be  locally 
defined  in  terms  of  the  measurable  parameters  of  temperature, 
salinitv,  speed  of  sound  and  surface  wave  height. 

(c)  the  mounting  of  the  hydrophone  be  stiff  enough  to 
preclude  significant  changes  of  sound  amplitude  through  the  sound 
beam. 

(d)  the  path  length  be  long  enough  to  observe  measurable 
amplitude  modulation. 

The  compromise  transducer  to  hydrophone  distance  was  set 

at  two  meters.  Checking  requirement  (a),  for  the  available  source 

2 

(F-27)  using  the  criterion  r  =  a  A,  and  assuming  an  average  sound 
speed  of  1500  meters  a  second; 

f  =  60  x  10^  Hz 

X  =  c/f  =  (  1 500/6 ) x  10-/+=  2.50  x  10-2  meters 
a  =  0.127  meters,  where  a  is  the  radius  of  the 
transducer  face 

15 


du: 
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then,  the  far  field  begins  at  (approximately) 

r  -  a2 A  =  0.01613/0.0250  =  0.645  meters 

Thus  the  far  field  begins  veil  within  the  two  meter  dis¬ 
tance  between  the  transducer  and  the  hydrophone. 

For  requirements  (b)  and  (d)  the  two  meter  distance  was 
selected  to  be  long  enough  to  yield  measurable  amplitude  modu¬ 
lation,  yet  short  enough  to  allow  oceanographic  measurements  on  a 
definable  volume  of  water.  The  data  from  the  runs  proved  this 
selection  to  be  justified. 

The  alignment  of  the  receiving  hydrophone  with  the  acoustic 
axis  of  propagation  was  an  important  consideration.  A  narrowly 
defined  beam  pattern  reduced  reflection  from  the  apparatus,  but  a 
radius  of  error  had  to  be  calculated  in  order  to  establish  the 


effect  of  minute  movement  of  the  receiver  (requirement  (c)). 


Pressure  amplitude  as  a  function  of  angle  was  calculated 
from  the  formula  below,  where  only  the  last  term  is  important; 


£(<Ail-4nr)  r  z  X,  ( not  if*.**' 


K  &  0 


For  one  per  cent  modulation,  an  angular  displacement  of 
0.45°  would  have  to  occur,  two  meters,  the  radius  of  movement 
would  be; 

9  =  0.45° 

sin  9  =  x/2  meters 

x  =  2(0.0075) 

x  =  0.0150  meters 

The  radius  of  movement  would  have  to  be  1.50  centimeters, 
more  movement  than  the  mounting  would  allow.  As  a  final  test  at 
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the  scene  of  the  experiment,  a  div^r  "twanged"  the  hydrophone 
supporting  wire  with  no  observable  result  in  amplitude  modulation. 

It  is  estimated  that  no  more  than  0.19Z  modulation  could  be  caused 
by  movement  of  the  hydrophone  by  ocean  forces. 

4 •  Design  of  the  Apparatus 

The  apparatus  (see  Figure  II-2)  as  designed  was  a  square 
of  height  six  feet  and  length  two  meters  from  the  face  of  the 
transducer  to  the  receiving  hydrophone. 

Thtee  sides  of  the  mount  were  of  two  inch  diameter  steel 
pipe  and  the  fourth  was  one-eighth  inch  steel  wire  stretched  by 
means  of  steel  springs  to  a  tension  of  one  hundred  and  fifty 
pounds . 

The  square  shape  offered  two  distinct  advantages,  the 
shortness  of  the  horizontal  members  minimized  vibratory  "tuning 
fork"  effects  in  conjunction  with  the  tension  in  the  wire;  it 
also  minimized  the  number  of  reflection  paths  from  source  to 
receiver  Acoustic  absorbent  rubber  (SoAB)  mounted  for  eighteen 
inches  along  the  upper  and  lower  members  also  reduced  possible 
reflections  to  the  point  that  reflections  were  down  38  db  relative 
to  the  signal  when  a  pulsed  60  kHz  signal  was  propagated. 

An  eleven  inch  steel  disc  was  welded  to  the  vertical  pipe 
member  directed  perpendicularly  to  the  sound  axis  as  a  mount  for 
the  transducer . 

The  completed  design  combined  two  necessary  characteristics, 
rigidity  and  minimum  self-generated  hydrodynamic  turbulence.  The 
end  of  the  apparatus  holding  the  transducer  was  necessarily 
massive  to  insure  rigidity  about  the  p.opagating  transducer. 
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ARATUS  IN  EXPERIMENTAL  CONFIGURATION 
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Hydrodynamic  turbulence  scattered  from  the  transducer  was  kept  to 
a  minimum  by  having  the  sound  source  end  of  the  entire  apparatus 
face  the  mean  current. 

The  wire  served  as  the  mounting  for  the  receiving  hydro¬ 
phone.  Since  the  hydrophone  was  quite  light  in  weight,  the  wire 
was  rigid  enough,  while  the  small  size  of  the  wire  precluded  any 
significant  hydrodynamic  turbulence. 

5 •  The  Transducer 

The  transducer  utilized  was  the  USRD  type  F-27.  The 
transmitting  array  consists  of  fifty-five  one-inch  diameter  lead 
metaniobate  discs  0.22  inches  thick,  each  of  which  is  cemented  to 
a  one-inch  diameter,  ^  inch  thick  tungsten  backing  plate  which 
acts  approximately  as  a  circular  piston. 

The  F-27  transducer  was  designed  for  wide  range  unidi¬ 
rectional  underwater  sound  projection  at  a  frequency  range  of  one 
kilohertz  to  forty  kilohertz.  The  USRD  recommends  a  limit  of 
forty  kilohertz  since  the  beam  pattern  was  considered  to  be  too 
sharp  at  higher  frequencies.  This  characteristic  was  advantageous 
to  the  experiment  at  sixty  kilohertz  and  results  were  indicative 
of  proper  transducer  selection,  with  the  one-half  power  points 
measured  at  three  degrees  off  the  sound  axis,  and  the  side  lobes 
at  sixty  degrees  down  forty  db  from  the  axial  sound  pressure  level. 

6 .  The  Receiving  Hydrophone 

The  experiment  required  a  small,  sensitive  hydrophone.  The 
hydrophone  used  was  the  Atlantic  Research  Corp.  type  LC-10  with  a 
diameter  of  0.635  centimeters.  For  the  sixty  kilohertz  carrier 


wave,  the  wavelength  is  2.50  cm,  making  the  wavelength  3.95  times 


the  diameter  of  the  hydrophone.  The  sensitivity  of  the  LC-10  is 
-110  dB  re  1  volt  per  microbar  at  60  kHz.  The  hydrophone  was 
mounted  vertically  along  the  wire,  with  the  central  part  of  the 
cylinder  perpendicular  to  the  sound  axis. 

B.  PROCESSING  ELECTRONICS 

1 .  Overall  System  Operation 

The  processing  electronics  were  designed  to  present  the 
waveform  of  the  amplitude  modulation  taking  place  on  the  acoustic 
carrier  wave  over  a  dynamic  voltage  range  compatible  with  that  of 
the  recording  equipment.  Generally,  it  was  desired  to  present  a 
ten  volt  peak-to-peak  signal  at  the  entry  leads  of  the  envelope 
detector,  and  then  process  the  waveform  to  view  only  the  amplitude 
modulation.  (See  Figure  II-3-) 

The  60  kHz  wave  was  originated  in  a  GR  1162A  frequency 
synthesizer  which  produced  a  one  volt  to  four  volt  peak-to-peak 
waveform  accurate  and  constant  to  within  one  millivolt  over  ex¬ 
tended  periods  of  operation. 

The  carrier  waveform  was  then  amplified  to  the  voltage 
necessary  to  produce  a  ten  volt  peak-to-peak  waveform  at  the 
envelope  detector  by  an  HP  467A  power  amplifier  set  at  a  fixed 
"times  ten"  setting.  It  was  found  to  be  advantageous  to  set  the 
power  amplifier  at  a  fixed  output  and  vary  the  output  of  the 
GP  1162a  frequency  synthesizer  in  order  to  achieve  the  desired 
voltage  value  at  the  envelope  detector. 

The  power -amp 1  if ied  signal  was  then  transmitted  through 
one  hundred  and  fifty  feet  of  waterproof  cable  to  the  F-27  trans¬ 
ducer  and  transmitted  to  the  LC-10  receiving  hydrophone. 
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ELECTONICS  PROCESSING  AND  CATA  GATHERING  SYSTEM  CONFIGURATION 


The  signal  as  received  by  the  LC-10  hydrophone  was  amp¬ 
lified  by  a  fixed  +30  db  (x  31*6)  gain  NUS  model  2010-030  voltage 
preamplifier  in  a  watertight  can  mounted  on  the  upper  horizontal 
member  of  the  mounting  apparatus.  The  preamplifier  required  a 
twelve  volt  DC  input,  supplied  by  a  Hewlett-Packard  721A  DC  power 
supply.  The  amplified  signal  was  transmitted  through  one  hundred 
and  fifty  feet  of  watertight  cable  to  the  initial  electronic  pro¬ 
cessing  equipment. 

The  initial  processing  of  the  modulated  waveform  consisted 
of  bandpass  filtering  the  signal  with  a  center  frequency  of  60  kHz 
through  a  pass  band  of  +_  three  hundred  Hertz,  with  -24  db  per 
octave  attenuation  on  either  side  of  the  pass  band.  The  unit 
utilized  was  the  Khron-Hite  model  3322  variable  filter.  In  order 
to  arrive  at  a  ten  volt  peak-to-peak  waveform  at  the  output  of 
the  filter,  a  +20  db  internal  amplification  feature  of  the  filter 
was  utilized.  The  output  at  this  point  consisted  of  a  ten  volt 
maximum  peak-to-peak  voltage  with  a  low  frequency  environmental 
amplitude  modulation  superimposed  upon  it,  varying  to  a  maximum  of 
sixteen  per  cent  of  the  peak  wave  form. 

The  above  waveform  was  processed  to  detect  the  modulation 
by  the  envelope  detector  (btc  Figure  II-4)  described  in  detail  in 
the  following  section.  The  salient  point  is  that  the  ten  volt 
signal  was  reduced  by  a  factor  of  0.6  in  the  audio  transformer  of 
the  envelope  detector. 

The  output  of  the  envelope  detector  was  the  varying  voltage 
induced  by  the  amplitude  modulation  superimposed  upon  a  small 
amount  of  DC  voltage. 
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The  envelope  detector  output  was  led  into  the  positive 
input  jack  of  a  PAR  Type  113  differential  amplifier.  The  dif¬ 


r 


ferential  amplifier  was  used  to  eliminate  the  IX  component  of  the 
amplitude  modulation  waveform.  This  was  accomplished  by  simultan¬ 
eously  placing  a  DC  voltage  into  the  negative  input  jack.  The  DC 
voltage  was  supplied  by  an  HP  721A  DC  power  supply  set  at  an 
output  of  four  volts  and  reduced  by  means  of  a  500  ohm  potentio¬ 
meter  to  that  (negative)  voltage  required  to  reduce  the  amplitude 
modulation  to  a  purely  alternating  signal. 

The  potentiometer  was  necessary  to  assure  a  constant  DC 
input  level.  The  power  supply/potentiometer  combination  was 
tested  over  a  continuous  two  hour  period  with  less  than  a  one 
millivolt  variation  during  the  entire  period.  The  potentiometer 
was  configured  with  an  outsized  tuning  wheel  to  allow  minute  (less 
than  one  millivolt)  adjustments  of  the  DC  compensating  voltage. 
(See  Figure  I I- 5 . ) 

The  compensated  amplitude  modulated  signal  was  amplified 
by  a  factor  of  ten  in  the  differential  amplifier  which  made  the 
varying  voltage  too  high  for  the  tape  recorder,  but  which  was 
subsequently  reduced  by  a  potentiometer.  The  ten-times  setting 
was  the  minimum  amplitude  setting  on  the  PAR  113* 

The  signal  was  subsequently  passed  through  another  Khron- 
Hite  3322  filter  set  in  a  four -pole  Butterworth  low  pass  configu¬ 
ration  with  a  high  cutoff  of  ten  Hertz  and  a  minimum  frequency  of 
0.001  Hertz.  The  attenuation  at  the  high  frequency  end  was  -48  db 
per  octave.  The  internal  amplification  feature  was  not  utilized. 
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The  lowpass -f il ter ed  amplitude  modulation  signal  was  re¬ 
duced  to  allow  dynamic  voltage  range  compatibility  with  the  tape 
recorder.  The  reduction  was  accomplished  by  the  use  of  a  250  ohm 
potentiometer  set  to  reduce  the  low-pass  filter  output  by  a  factor 
of  ten.  This  reduction  had  the  desirable  result  of  restoring  the 
amplitude  modulation  t.o  its  value  relative  to  the  six  volt  output 
of  the  audio  transformer  in  the  envelope  detector. 

The  signal  was  recorded  on  a  Sangamo  model  3500  fourteen 
track  tape  recorder  simultaneously  with  the  oceanographic  mea¬ 
surements  during  each  run.  The  dynamic  range  of  the  tape  recorder 
in  the  FM  mode  utilized  was  0.01  to  2.50  volts,  with  all  the  amp¬ 
litude  modulation  observed  falling  within  this  range. 

The  electronic  processing  design  proved  to  be  quite 
satisfactory  except  for  details  in  various  components  which  were 
amply  compensated  for  during  the  course  of  the  experiment. 

2 .  The  Envelope  Detector 

The  envelope  detector  was  adapted  from  an  ocean  wave  height 
detector  and  proved  to  be  satisfactory.  The  detector  is  essentially 
a  double  rectification  device  with  an  RC  filter  and  a  capactive  lead 
to  ground  potential  to  discnarge  DC  voltages  built  up  in  the  sili¬ 
con  diode  ( 1N540)  rectifying  elements.  (See  Figure  II-4.) 

The  time  constant  for  the  RC  filter  at  a  carrier  fre¬ 
quency  of  60  kHz  is  fifty  microseconds  and  with  double  rectifi¬ 
cation  there  was  no  appreciable  decay  between  peaks. 

The  envelope  detector  proved  sensitive  to  acoustic  modu¬ 
lation  of  as  little  as  one  twentieth  of  one  per  cent  in  the 
laboratory  tanks. 

2  5 
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Possible  improvements  to  the  detector  would  be  the  sub¬ 
stitution  of  larger  capacitors  in  the  filtering  element  to  reduce 
the  substantial  (1.3  volt)  DC  component  in  the  output. 
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III. 


EXPERIMENTAL  PROCEDURE 


A.  GENERAL 

The  amplitude  modulation  runs  were  performed  on  the  twenty- 
first  of  October  1971  at  the  NUC  Oceanographic  Research  Tower 
at  San  Diego.  The  runs  were  twenty  minutes  in  length,  to  allow 
enough  time  for  major  low  frequency  ocean  phenomena  to  be  ade¬ 
quately  described  and  measured. 

The  ocean  sensors  were  mounted  around  the  acoustic  sensors 
in  o  • er  to  define  the  volume  of  ocean  over  which  the  amplitude 
modulation  took  place.  (See  Figure  1 1 1  —  1 . )  The  sensors  were 
selected  to  measure  temporal  variations  of  temperature,  salinity, 
turbulent  veloci ty,  wave  height  and  the  speed  of  sound  during  the 
period  of  acoustic  measurement.  During  each  run,  the  oceanogra¬ 
phers  present  also  measured  and  observed  weather  and  sea  surface 
conditions  as  background  for  the  more  precise  information  pre¬ 
sented  by  the  ocean  sensors.  The  ocean  sensors  used  are  listed 
below  with  respect  to  the  parameters  they  measured. 

1 .  Temperatur  e 

a.  Three  thermistors  located  about  the  apparatus  obtained 
temporal  variation  of  temperature  in  the  volume  of  water  under 
obser  vat  ion . 

b.  One  Bisset -Ber man  temperature  sensor  measured  absolute 
values  of  temperature  near  the  F-27  transducer  over  a  range  of 
14°C  to  19°C. 
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Figure  III- 
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One  Bisset-Berman  salinometer  located  near  the  F-27 
transducer  measured  temporal  variation  in  salinity  over  a  range 
of  37.5  to  39.5  parts  per  thousand. 

3 •  Particle  Velocity 

An  electromagnetic  flowmeter  measured  one  component,  (u) 
or  (v),  of  horizontal  velocity  and  the  vertical  component  (w).  The 
flowmeter  was  located  about  one-half  of  the  distance  between  the 
acoustic  sensors  1.13  meters  downstream  from  the  acoustic  axis. 

4 .  Speed  of  lound 

The  sound  speed  was  measured  as  an  absolute  value  with 
temporal  variations  by  a  Ramsey  Corp.  Mark  I  SVTD  probe  located 
near  the  base  of  the  acoustic  apparatus. 

5.  Wave  Height 

The  wave  height  was  measured  by  two  sensors,  a  pressure 
wave  height  indicator  and  a  Baylor  gauge. 

As  noted  above,  not  all  the  ocean  sensors  were  operating 
during  all  runs.  Five  runs  were  made,  with  runs  two  through  four 
considered  suitable  for  analysis.  The  general  scheme  was  to  per¬ 
form  two  runs  near  the  surface,  the  first  run  with  the  electro¬ 
magnetic  flowmeter  measuring  in  the  u  and  v  (horizontal) 
directions,  and  the  second  with  the  flowmeter  measuring  in  the  u 
and  w  (or  vertical)  directions.  Three  runs  followed  at  deeper 
depths  with  the  flowmeter  in  the  u-w  configuration.  The  first 
run  was  aborted  when  the  flowmeter  malfunctioned,  but  runs  two 
through  five  proved  to  be  successful.  A  listing  of  the  runs  taken 
is  presented  below;  all  times  are  local  and  the  date  is  21  October 
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1971.  Mean  water  level  was  6.7  meters  below  the  concrete  deck 
of  the  NUC  tower.  The  depths  are  in  meters  below  mean  water. 

Run  1 . 

depth:  4.3  meters, 

times:  start:  1419  stop:  1439 

comment:  The  data  from  this  run  were  not  analyzed  due  to 
a  malfunction  of  the  electromagnetic  flowmeter. 

Run  2. 

depth:  4-3  meters , 

times:  start:  153°  stop:  1550 

sensors;  turbulence:  1.  u  (horizontal)  velocity 

2.  w  (vertical)  velocity 

3.  sound  velocity 

4.  Both  wave  heights 
temperature:  1.  thermistor  T^ 

2.  Bissett -Berman  temperature 
salinity:  1.  Bissett-Ber man  salinometer 
Run  3  • 

depth:  9.3  meters, 

times:  start:  l6l6stop:  1636 

sensors:  turbulence:  1.  u  velocity 

2.  w  velocity 

3.  sound  velocity 

4.  both  wave  heights 

temperatures:  1.  thermistors  T^,T2,  and  T^ 

2.  Bissett-Berman  temperature 
salinity:  Bissett-Berman  salinometer 
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Run  4- 


depth:  14.3  meters, 

times:  start:  1 6 48  stop:  1708 

sensors:  all  sensors  operating,  see  sensor  information 
for  Run  3- 
Run  5 . 

depth:  7-3  meters, 

times:  start:  1728  stop:  1748 

sensors:  all  sensors  operating,  see  sensor  information 
for  Run  3- 

B.  CONFIGURATION  OF  THE  OCEAN  SENSORS  IN  RELATION  TO  THE  ACOUSTIC 
APPARATUS 

Two  factors  determined  the  configuration  of  the  ocean  sensors 
relative  to  the  acoustic  apparatus.  First,  the  volume  of  water 
around  the  acoustic  axis  had  to  be  described  adequately  for 
temporal  variations,  and  second,  the  ocean  sensors  had  to  be 
placed  such  that  there  was  little  or  no  interference  from  down¬ 
stream  turbulence  from  the  sensors  themselves. 

Temperature  variations  and  turbulent  velocities  were  considered 
to  be  the  most  important  parameters  affecting  the  amplitude  modu¬ 
lation.  Temperature  variations  were  recorded  by  the  three  ther¬ 
mistors  and  the  Bi ssett -Ber man  temperature  probe.  The  Bissett- 
Berman  probe  was  placed  centrally  between  the  acoustic  apparatus 
and  the  mount  for  the  flowmeter  due  to  its  size.  Two  of  the  ther¬ 
mistors  were  placed  above  and  below  the  receiving  hydrophone  in  the 
vertical  plane  of  the  acoustic  axis,  while  the  third  was  placed 
along  the  horizontal  plane  of  the  axis  midway  between  the  trans¬ 
ducer  and  the  hydrophone.  (See  Figure  III-2.) 
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Figure  III- 


The  electromagnetic  flowmeter  was  placed  so  that  it  measured 
the  turbulence  occurring  113  centimeters  from  the  center  of  the 
acoustic  axis . 

All  other  sensors  were  integral  parts  of  large  instruments 
and  were  placed  between  the  acoustic  apparatus  and  the  mount  for 
the  flowmeter  to  minimize  the  effect  of  their  wake  turbulence. 

Figure  1 1 1  —  3  shows  the  distances  between  various  sopors.  The 
mean  current  was  as  shown  in  Figure  III-l,  with  the 
apparatus  up-stream  of  all  other  equipment. 

C.  DETERMINATION  OF  RUN  DEPTHS 

The  mean  water  depth  at  the  NUC  tower  was  16.4  meters  (54  feet) 
during  the  time  of  the  runs.  Four  depths  were  selected  in  order  to 
observe  near-surface  to  near-bottom  amplitude  modulation  effects. 
The  first  two  runs  were  at  an  estimated  water  depth  of  4*3  meters, 
the  third  at  an  estimated  9.3  meters,  the  fourth  (and  deepest)  at 
14.3  meters,  and  the  last  run  was  at  7.3  meters. 

These  depths  were  chosen  in  order  to  cover  the  entire  water 
column  and  to  minimize  interference  from  the  underwater  beam 
structure  of  the  tower.  The  shallowest  depth  was  far  enough  from 
the  surface  to  avoid  undue  stresses  on  the  apparatus,  while  the 
deepest  depth  was  far  enough  away  from  the  bottom  to  avoid  inad- 
vertant  contact. 

D.  THE  NUC  TOWER 

The  NUC  tower  is  located  approximately  one  mile  off  Mission 
Beach,  California  in  16  meters  of  water.  The  apparatus  as  de¬ 
picted  in  Figure  III-2  was  mounted  on  the  west  side  of  the  tower. 
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Figure  III- 


The  electronics  were  installed  on  the  main  deck,  while  the  appa¬ 
ratus  was  handled  from  the  second  deck  immediately  below.  (See 
Figure  III -4  -  ) 

Environmental  factors  were  observed  at  the  tower  while  the 
runs  were  in  progress.  The  factors  noted  were: 

1.  swell:  from  the  northwest,  two  feet  from  crest  to 
trough,  with  an  estimated  period  of  four  to  five  seconds. 

2.  wind:  310°,  300°,  320°  from  9  knots  to  7.5  knots  during 
the  course  of  the  experiment 

3.  wind  waves  six  inches  from  crest  to  trougn 

4.  some  bubbles  visible  on  the  surface 

E.  UNUSUAL  EFFECTS 

Unusual  effects  were  those  effects  arising  from  considerations 
other  than  physical  parameters  and  naturally  occurring  bubble 
populations.  All  the  unusual  effects  noted  were  due  to  animal  and 
plant  life  in  and  around  the  apparatus  during  each  run. 

The  area  around  the  tower  supported  a  large  population  of 
small  fish  ranging  up  to  tight  inches  in  length.  It  was  noted 
that  fish  were  attracted  to  the  apparatus  when  it  was  lowered 
initially,  but  resumed  their  normal  habits  after  a  short  period 
of  ti.:ie.  Fish  could  not  affect  the  wire  mounted  hydrophone,  nor 
were  any  effects  noted,  since  fish  swimming  through  the  sound  beam 
would  affect  the  modulation  for  a  minimal  amount  of  time. 

Plant  life  could  affect  the  equipment  by  drifting  into  it. 
There  was  a  small  amount  of  kelp  in  the  area  at  all  depths,  but 
no  acoustic  effects  noted  during  the  runs. 
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Boat  traffic  around  the  tower  had  some  effect,  and  could  have 
been  sources  of  turbulent  velocity  and  bubbles  which  could  noti- 
cibly  affect  the  amplitude  modulation  of  the  acoustic  signal.  This 
effect  was  quite  apparent  during  the  Run  5,  when  there  was  a  boat 
idling  in  the  area. 

F.  ORGANIZATION  OF  INDIVIDUAL  RUNS 

All  information  from  each  sensor  was  recorded  on  Sangamo 
Model  3500  tape  recorder.  Eetween  individual  runs,  all  sensors 
were  monitored  for  proper  output.  After  a  depth  change,  a  few 
minutes  were  allowed  for  all  instruments  to  be  checked  for  proper 
operation.  When  all  instruments  were  checked  and  operating,  the 
tape  recorder  was  started.  The  simultaneous  start  and  stop  times 
controlled  by  the  tape  recorder  insured  an  accurate  run  period 
for  all  sensors. 

The  recorder  channel  allocation  was  as  follows: 

Channel  1.:  Ramsey  prcbe  temperature 
Channel  2.:  Ramsey  probe  velocimeter 
Channel  3> :  Bissett -Ber man  salinometer 
Channel  4.:  Turbulent  velocity  (v) 

Channel  5.:  Turbulent  velocity  (w) 

Channel  6.:  Baylor  wave  height  probe 
Channel  7-:  Amplitude  modulation 
Channel  8.:  Pressure  wave  height  probe 
Channel  9.:  Thermistor  T^ 

Channel  10.  :  Bissett-Berman  temperature 
Channel  11.:  Thermistor  T 2 
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Channel  12.:  Thermistor 
Channel  13*:  Phase  modulation 
Channel  14  - :  Voice 


IV.  ANALYSIS  OF  ACOUSTIC  RESULTS 


A.  INITIAL  DATA  REDUCTION 

1 .  General 

The  recorded  a  nr.  log  data  consisted  of  four  twenty  minute 
runs  at  varying  depths.  The  dynamic  range  of  the  recorded  voltage 
was  0.01  to  1.10  volts.  The  data  reduction  scheme  involved  three 
steps;  transfer  of  the  analog  data  from  magnetic  tape  to  a  Brush 
Mark  200  strip-chart  recorder,  conversion  of  the  printed  analog 
data  to  digitized  data  and  computer  analysis  of  the  digitized  data. 

2 .  Magnetic  Tape  to  Strip-Chart  Record 

The  tape  recorded  data  was  originally  collected  at  a  tape 
speed  of  1  7/A  inches  per  second.  When  replayed  for  the  Brush 
Mark  200  eight  channel  -strip-chart  recorder,  the  data  was  played 
at  7h  inches  per  second,  making  the  strip-chart  data  for  each  run 
one  fourth  as  long  as  the  real-time  runs.  This  had  the  effect  of 
compressing  a  twenty  minute  run  to  five  minutes  of  strip-chart 
readout  and  increasing  the  apparent  freouencies  of  modulation  by  a 
factor  of  four.  At  a  strip-chart  recorder  speed  of  five  milli¬ 
meters  per  second,  each  printed  run  was  presented  on  1.5  meters 
of  chart. 

The  strip-chart  recorder  was  set  to  present  the  data  at 
50  millivolts  per  millimeter,  or  1.27  volts  per  inch. 

The  other  seven  strip-chart  channels  were  used  to  visually 
compare  some  of  the  ocean  sensor  outputs  to  the  output  from  the 
amplitude  modulation  data.  The  other  seven  strip-chart  channels 
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were  allocated  as  follows: 

Channel  1.:  Salinometer 
Channel  2.:  v  velocity 
Channel  3*:  w  velocity 
Channel  b>-  Amplitude  modulation 
Channel  5*:  Baylor  wave  height 
Channel  6.:  Thermistor  T^ 

Channel  7.:  Bissett-Berman  temperature 
Channel  8.:  Velocimeter 

It  was  then  possible  to  make  preliminary  visual  comparison 
of  the  waveforms  of  the  recorded  outputs  in  order  to  observe  any 
similarities  in  the  waveforms.  {See  Figures  TV-1  through  IV-4. ) 

3 .  Strip-Chart  Record  to  Digital  Tape 

The  strip-chart  records  of  each  of  the  four  runs  to  be 
analyzed  were  taken  to  the  Fleet  Numerical  Weather  Central  facility 
at  Point  Pinos,  California  to  be  recorded  as  digitized  data  on 
magnetic  tape.  The  data  were  digitized  directly  from  the  strip 
chart  recording  on  the  facility's  tracing  digitizer.  The  digi¬ 
tizer  recorded  data  points  in  increments  of  hundredths  of  an  inch 
in  the  y-direction  for  each  one  one-hundredth  of  an  inch  of  'ravel 
in  the  x-direction,  and  transferred  the  x  ana  y  values  to  a  nine- 
track  magnetic  tape. 

The  digitizing  was  initiated  at  the  start  of  each  run, 
visible  on  the  strip-chart  as  a  vertical  line  produced  when  the 
tape  recorder  run  was  commenced.  The  zero  amplitude  value  in  the 
y-direction  and  the  zero  time  value  in  the  x-direction  were  taken 
at  this  line.  There  were  data  excursions  in  the  positive  and 
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EIGHT  CHANNEL  BRUSH  RECORDER  READOUT  FOR  RUN  2,  DEPTH: 
4.3  meters,  top  to  bottom:  sa 1 inometer ,  v-turbulence , 
w-tur bulence ,  amplitude  modulation,  wave  height,  ther¬ 
mistor  3,  Bissett -Berman  temperature  and  velocimeter 
readouts  . 


Figure  I V - 1 


EIGHT  CHANNEL  BRUSH  RECORDER  READOUT  FOR  RUN  5  DEPTH: 
7.3  meters.  For  channel  identification,  see  Fig.  IV-1 

Figure  IV-2 


EIGHT  CHANNEL  BRUSH  RECORDER  READOUT  FOR  RUN  3  DEPTH: 
9.3  meters.  For  channel  identification,  see  Fig.  IV-1 


Figure  IV-3 


EIGHT  CHANNEL  BRUSH  RECORDER  READOUT  FOR  RUN  4  DEPTH: 
14.3  meters.  For  channel  identification  see  Fig.  IV-1 


Figure  IV-4 


negative  directions,  measured  in  one-hundredths.  The  conversion 
factors  were  0.0127  volts  per  data  point  in  the  y-direction,  and 
0.01265  seconds  (real  time)  per  data  point  in  the  x-direction. 

4 •  Digitized  Data  from  Magnetic  Tape  to  IBM  Data  Cards 

The  conversion  of  the  digitized  data  from  magnetic  tape 
to  data  cards  for  analysis  on  the  IBM  360  digital  computer  was 
accomplished  in  the  Fleet  Numerical  Weather  Central's  CDC  5000 
computer . 

The  procedure  utilized  Program  I,  and  converted  at  least 
5600  data  points  for  each  run,  corresponding  to  nineteen  minutes 
of  real  data  collection  time.  This  decrease  in  data  time  was  due 
to  an  uncertainty  of  stopping  points  on  the  strip-chart  readouts, 
and  allowed  an  identical  number  of  data  points  for  each  run. 

B.  ANALYSIS  OF  ACOUSTIC  DATA 

The  data  from  each  run  were  analyzed  on  the  Naval  Postgraduate 
School's  IBM  360  digital  computer. 

The  analyses  were  performed  on  time-varying  amounts  of 
fractional  amplitude  modulation.  The  raw  data  as  presented  on 
data  cards  was  in  varying  values  expressed  as  positive  and 
negative  hundredths  of  an  inch.  The  conversion  was  common  to  all 
programs  and  is  described  below: 

a.  Each  run  was  scanned  for  the  maximum  positive  value. 

b.  The  maximum  positive  value  was  subtracted  from  all  other 
values  and  converted  to  0.0.  This  level  corresponded  to  the 
unmodulated  amplitude.  All  data  values  were  now  negative  values, 
all  less  than  zero  except  for  the  maximum  value. 
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c.  Each  value  was  multiplied  by  a  conversion  factor  consisting 


of  (-1.27)  which  converted  the  data  points  to  positive  voltage 
values,  and  (1/6)  which  made  each  point  a  fractional  value  of  the 
maximum  amplitude  of  six  volts.  The  resultant  data  points  could 
now  be  expressed  as  fractional  modulation: 


F 


(unmodulated  voltage)  -  (modulated  voltage) 
(unmodulated  voltage) 


The  converted  raw  data  was  then  analyzed  for  the  mean, 
variance,  temporal  correlation  and  power  spectral  density  of  each 
run  in  Program  II.  Program  III  computed  a  histogram  of  the  amount 
of  fractional  modulation  versus  the  number  of  occurrences  of  each 
value.  An  independent  calculation  of  the  mean,  mea  square  and 
variance  was  made  in  Program  IV.  (See  Figures  IV-5  through  IV-7.  ) 
The  results  of  the  various  analyses  are  presented  in  the  following 
section. 
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DATA  ANALYSIS  SCHEME:  PROGRAM  II 


AT  THIS  POINT,  ALL  VALUES  MAY  BE  EXPRESSED  AS: 

.  _  (UNMODULATED  VOLTAGE)  -  (MODULATED  VOLTAGE) 
(UNMODULATED  VOLTAGE) 
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Figure  IV-5 


DATA  ANALYSIS  SCHEME:  PROGRAM  III 


RAW  DATA  ON  IBM  CARDS 


CONVERSION  OF  ALL  DATA  VALUES 
TO  FRACTIONAL  MODULATION  VALUES 


ESTABLISH  "BINS"  FROM  0.0  TO  0.170 
IN  STEPS  OF  0.001,  YIELDING  171  "BINS" 
(0.160  IS  THE  MAXIMUM  VALUE  OF 
FRACTIONAL  MODULATION) 


± _ 

SCAN  EACH  VALUE  OF  FRACTIONAL  MODULATION 
AND  PLACE  IN  THE  APPROPRIATE  BIN 


_ ± 

COUNT  THE  NUMBER  OF  OCCURENCES  IN  EACH  BEN 


Figure  IV-6 


Figure  IV-7 
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V.  ACOUSTIC  RESULTS 


A  THE  STRIP-CHART  READOUTS 

Figures  V-l  through  V-4  are  representative  168-second  portions 
of  each  twenty-minute  run.  The  variation  in  voltage  from  peak-to- 
peak  is  approximately  one  volt,  and  the  apparent  frequencies  of 
variation  are  four  times  the  actual  frequencies  occurring  during 
each  run. 

From  visual  inspection  of  the  strip-charts,  it  may  be  seen 
that  the  modulations  occurring  during  the  two  shallow  runs  (4-3 
and  7.3  meters)  were  more  active  in  frequency  than  during  the  two 
deeper  runs  (9-3  and  l4 . 3  m.).  It  must  be  noted  that  the  7.3  meter 
run  was  conducted  after  the  other  runs,  and  that  during  most  of  the 
period  a  fifty-foot  motor  launch  was  in  the  area  of  the  tower  with 
with  propellors  turning.  The  launch  served  as  a  bubble  injection 
device  as  well  as  a  turbulence  generator.  The  ;un  was  not  invali¬ 
dated  because  the  modulation  effect  of  bubbles  artificially  gen¬ 
erated  during  this  run  could  be  compared  with  that  of  the  ambient 
medium  of  the  previous  runs. 


B.  MEANS  AND  VARIANCES 

The  mean  and  variance  of  the  amplitude  modulation  for  each  run 
was  determined  by  Computer  Run  IV.  The  formulation  was  based  on 
the  fractional  modulation  F,  and  was  calculated  as  shown: 


F  = 


VMAX  -  VMOD 


4K 


VMAX 


BRUSH  RECORDER  READOUT  OF  AMPLITUDE  MODULATION  RUN  TWO  DEPTH: 
4-3  meters,  vertical  lines  indicate  four  seconds  in  real  time 
horizontal  lines  indicate  0.25  volts. 
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BRUSH  RECORDER  READOUT  OF  AMPLITUDE  MODULATION  RUN  5  DEPTH: 
7.3  meters,  see  caption  for  Fig.  V-J . 

Figure  V-2 

•V* 


BRUSH  RECORDER  READOUT  OF  AMPLITUDE  MODULATION  RUN  3  DEPTH: 
9.3  meters,  see  caption  for  Fig.  V-l. 


BRUSH  RECORDER  READOUT  FOR  AMPLITUDE  MODULATION  RUN  4  DEPTH 
14.3  meters,  see  caption  for  Fig.  V-l. 


Figure  V-4 


which  resulted  in  the  time  series  F(t);  then, 

< F U)>  »  2  F; 


where  n  =  5600 


-- 


The  results  obtained  from  each  run  were  compared  with  the 
results  of  Computer  Run  II,  which  determined  the  mean,  temporal 
correlation  and  the  power  spectral  density  for  each  run.  Effecti¬ 
vely,  for  Program  II  the  mean,  (F(t)),  was  determined  as  in 

Program  IV  so  that  comparison  was  super flous, 

2 

The  variance,  O’  ,  can  be  formulated  in  three  ways;  first  as 
shown  above,  secondly  as  the  zero-time  temporal  correlation  value, 

t  T/* 

6  =  c<« j  -iuu.jp.  J  F({)F/t*rUi 
I*.  T-f  <*>  T 

,eo 

and  thirdly  as  the  value  of  the  area  under  the  curve  for  the  power 
spectral  density: 


<£1-- 


It  should  be  noticed  that  at  all  depths,  the  three  values  of 

the  variance  are  quite  close.  The  slightly  different  values  of 
2 

<7  from  the  calculation  for  C(0)  and  from  the  simple  sorting  cal¬ 
culation  (Program  IV)  are  solely  due  to  round  off  error  for  the 
latter  calculation. 
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The  value  of  O'  derived  from  the  integral  of  the  power  spectral 
density  has  the  largest  numerical  value.  This  is  the  r  esult  of 
limiting  the  higher  frequencies  from  the  calculation  of  the  Fourier 
pass  band,  resulting  in  a  rounded  "window".  This  rounding  off 
resulted  in  a  very  slightly  la  ger  area  under  the  power  spectral 
density  curve  than  would  otherwise  occur. 

The  results  for  each  run  are  shown  below: 

Run  2:  4.3  meters 

Program  IV:  < F>  =  0.0526 
<F2)  =  0.0041 
cr2  =  0.0013 

Program  II:<F>  =  0.05261 
ff2  =  c(0)  =  c. 00131 
a2  =  J  e ( f ) dF  =  0.0014469 

Run  5:  7*3  meters 

Program  IV:  < F>  =  0.0460 
<F2)  =  0.0031 

a2  =  0.0010 
Program  II:  ( F)  =  0.04604 
02  =  c(0)  =  0.00096 
O  = j  E(F)dF  =  0.0009756 
Run  3:  9.3  meters 

Program  IV:  < F>  =  O.O3O3 
<F2)  =  0.0016 

a2  =  0.0007 

Program  II:  < F>  =  0.03026 
a2  =  C(0)  =  0.00067 
a2  =  j  E( F) dF  =  0.0007378 
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Run  4:  i4  •  3  meters 

Program  IV:  <F>  =  O.0470 

<f2>  =  0.0040 
a2  =  0.0018 

Program  II:  (f)  =  0.0472 

a2  =  c(0)  =  0.00173 

O'2  =  i  E (  F ) dF  =  0.0019268 

The  probability  density  of  fractional  modulation  values  for 
each  run  is  shown  on  the  histograms.  (Figures  V-5  through  V-7.  ) 
Although  the  mean  values  calculated  above  do  not  differ  in  any 
orderly  manner  with  depths,  the  histograms  (as  well  as  Figures 
V-1,2,3  and  4)  show  that  for  the  shallow  depths  there  are  greater 
numbers  of  short  duration,  high  modulation  values,  whereas  for  the 
greater  depths,  the  distribution  of  fractional  modulation  centers 
heavily  about  the  lower  values.  The  bimodal  P.D.F.  at  the  greatest 
depth  ( 14.3  meters)  shov/s  that  onset  of  infrequent  occurrences  of 
very  high  modulation  values.  Iris  is  the  statistical  representation 
of  the  type  of  modulation  observed  on  Figure  V-4,  where  the  modu¬ 
lation  assumes  a  high  value  and  varies  about  that  value  for 
relatively  long  periods.  The  rla’a  thereby  suggest  that  as  depth 
is  increased  the  magnitude  of  the  medium  irihomogenieties  increase 
and  the  size  (temporal  extent)  of  the  inhomogeneous  volumes  becomes 
greater  as  well. 

C.  TEMPORAL  CORRELATION 

The  temporal  correlations  (time  autocorrelations)  for  each  run 
are  shown  on  Figures  V-8  through  V-ll.  The  times  to  50% 
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DISTRIBUTION  OF  FRACTIONAL  MODULATION 


DISTRIBUTION  OF  FRACTIONAL  MODULATION 
Figure  V-7 
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DEPTH:  4.3  meters 

Figure  V-8 
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NORMALIZED  TEMPORAL  CORRELATION  OF  RUN  5 
DEPTH:  7.3  meters 


Figure  V-9 


Fractional  Correlation 


NORMALIZED  TEMPORAL  CORRELATION  OF  PUN  3 
DEPTH:  9.3  meters 
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Figure  V-10 


Fractional  Correlation 


NORMALIZED  TEMPORAL  CORRELATION  OF  RUN  4 
DEPTH:  14-3  meters 

Figure  V-ll 
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correlation  and  1/e  correlation  a^e  shown  below: 


Time  to  1/e: 


Depth  (meters) 

Time  (se 

4.3 

2.33 

7.3 

1.43 

9.3 

2.23 

14.3 

6.9 

4.3 

3-25 

7.3 

1 .65 

9.3 

2.73 

14.3 

11.75 

It  can  be  seen  from  the  figures  that  there  is  a  rapid  initial 
decorrelation  down  to  low  as  10^  for  depths  of  4-3  and  9.3  meters, 
and  a  slower  decorrelation  for  the  14.3  meter  depth.  All  of  these 
runs  show  a  fluctuation  between  steadily  decreasing  minima  and 
maxima  with  a  delayed  zero  crossing.  Whitmarsh,  Skudrzyk  and 
Urick  (1957)  hypothesized  that  such  fluctuation  was  an  indication 
of  the  size  of  the  inhomogeneous  "patches"  around  the  receiver. 

If  Whitmarsh,  Skudrzyk  and  Urick’s  hypothesis  is  correct,  then  the 
apparent  "patch"  sizes  show  an  increase  with  depth.  If  this  is  the 
case,  then  small-scale  turbulence  may  tend  to  dominate  at  all  depths 
near  the  surface. 

The  temporal  correlation  of  the  7.3  meter  depth  is  of  interest 
since  it  is  known  that  an  artificially  derived  bubble  and  turbulence 
field  was  introduced  into  the  volume  under  consideration.  The  7.3 
meter  depth  shows  the  most  rapid  decorrelation  to  zero  of  all  runs, 
and  twice  as  frequent  fluctuations  about  zero  as  for  the  undisturbed 
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water  at  other  depths.  This  suggests  that  the  ambient  "patches" 
of  water  were  chopped  into  smaller  patches  by  the  boat's  propellor. 

Another  hypothesis  might  be  that  the  inhomogeneous  volumes 
are  more  or  less  "frozen"  and  that  the  spatial  displacement  of 
these  volumes  about  the  acoustic  axis  may  cause  the  appearance 
of  temporal  variation.  Corroboration  of  these  hypothetical 
explanations  awaits  the  complete  analysis  of  the  data  gathered  by 
the  oceanographic  sensors. 

D.  POWER  SPECTRAL  DENSITY 

The  power  spectral  densities  (PSD)  of  each  run  are  plotted 
against  frequency  on  Figures  V-12  through  V-15.  The  power  spectial 
densities  show  the  highest  values  at  those  frequencies  (0.08  to 
0.12  Hz)  which  appear  to  correspond  with  the  observed  ocean  swell 
frequencies  prevailing  during  each  run. 

The  PSD's  show  a  diminishing  value  both  with  increasing  fre¬ 
quency  and  with  increasing  depth.  The  two  shallowest  depths 
(4.3  and  7.3  meters)  show  a  similarity  in  the  values  of  PSD  from 
0  to  1.0  Hz,  except  for  a  high  "DC"  value  in  the  plot  for  4.3 

meters,  resulting  in  a  higher  contribution  to  the  integrated  value 

2 

of  PSD  (O  )  at  that  depth.  A  similar  likeness  is  apparent  between 
the  two  greater  depths  (9.3  and  l4 . 3  meters).  For  frequencies 
greater  than  0.1  Hz,  the  PSD  at  the  two  shallower  depths  appear 
to  be  consistently  four  to  five  db  greater  than  for  the  deeper 
depths . 

The  frequencies  of  relative  maxima  along  each  of  the  PSD  curves 
show  similarities  at  each  end  of  the  spectrum.  At  the  low  fre¬ 
quency  end,  the  ocean  swell  frequencies  are  more  pronounced  at 
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POWER  SPECTRAL  DENSITY  OF  RUN  2  DEPTH:  4.3  meters 


POWER  SPECTRAL  DENSITY  OF  RUN  5  DEPTH:  7-3  meters 


POWER  SPECTRAL  DENSITY  OF  RUN  4  DEPTH:  l4 . 3  meters 


the  shallower  depths  of  4-3  and  7.3  meters.  These  maxima  are 
higher  for  the  7-3  meter  run  than  for  the  4-3  meter  run,  which 
may  again  shov  the  effect  of  the  boat  in  the  area  injecting 
bubbles  and  turbulence  into  the  medium.  At  9.3  and  14-3  meters, 
these  low  frequency  maxima  tend  to  dissappear  into  the  overall 
curve,  showing  a  decrease  of  the  component  of  the  amplitude  modu¬ 
lation  due  primarily  to  the  turbulent  effect  of  surface  waves. 

The  amplitude  modulation  at  the  higher  frequencies  (0.2  to 
1.0  Hz)  could  be  due  to  spatial  distributions  of  temperature, 
salinity  or  bubbles  which  are  seen  as  temporal  frequencies  due  to 
movement  across  the  acoustic  axis. 

E.  CONCLUSIONS 

The  conclusions  presented  here  await  the  additional  infor¬ 
mation  that  the  data  from  the  ocean  sensors  will  provide.  This 
information  will  be  available  in  the  theses  of  LCDR  Duchock  and 
LT's  Seymour  and  Bordy,  to  be  presented  at  the  Naval  Postgraduate 
School  in  March,  1972.  The  present  conclusions  are  based  upon 
the  acoustic  data  alone  and  are  listed  below: 

1.  Amplitude  modulation  near  the  ocean  surface  ranged  in 
value  from  zero  to  sixteen  per  cent.  The  maximum  amplitudes  were 
not  noticably  restricted  to  any  one  depth,  bat  as  depth  increased, 
the  distribution  of  modulation  values  tended  to  center  about  the 
lower  values. 

2.  The  temperal  correlation  curves  for  the  4*3>  9.3,  and  l4 . 3 
meter  runs  indicate  that  the  size  of  the  volumes  of  inhomogen iet ies 
("cells")  increases  as  the  depth  is  increased.  The  correlation 
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curve  for  the  7-3  meter  run  shows  a  significant  decrease  in  the 
size  of  the  volumes,  a  direct  result  of  having  a  motor  launch  in 
the  area  injecting  turbulence  and  bubbles  into  the  volume  of  water 
under  investigation.  Thus  it  may  be  seen  that  in  the  ambient 
ocean,  cells  of  microstructure  are  larger  at  greater  depth.  Near 
surface  cells  would  tend  to  be  reduced  in  size  by  the  action  of 
turbulence  arising  from  the  effect  of  surface  waves. 

3.  A  comparison  of  the  variances  of  the  two  shallow  runs, 
at  4-3  and  7.3  meters,  shows  that  there  is  a  small  difference 
between  the  two.  (CT2(4-3)  =  0.0013  and  CT2  ( 7  •  3 )  =  0.0010.)  Com¬ 
paring  the  temporal  correlations  for  the  two  runs  it  maybe  seen 
that  the  decorrelation  time  for  the  7-3  meter  run  is  one-half  that 
of  the  4*3  meter  run,  with  more  numerous  fluctuations  about  the 
zero  axis.  Such  data  suggests  a  reduction  in  the  size  of  the  cells 
of  inhoraogeniety . 

Mintzer's  formula  for  the  variance  of  acoustic  amplitude 
modulation,  derived  for  an  exponential  correlation  function  (see 


Appendix  I) 


shows  the  dependence  of  variance  on  the  wave  number,  k,  range,  r, 
the  RMS  variation  of  the  index  of  refraction,  0L ,  and  the  size  of 

the  inhomogen iet ies ,  a.  In  this  experiment,  wave  number  and  range 

2  2 

were  constant;  since  neither  a  nor  (CV)  changed  appreciably,  it 
can  be  deduced  that  there  was  a  simultaneous  increase  in  &  (with 
the  reduction  of  a)  for  the  7*3  meter  run. 

4.  The  parameters  which  could  cause  a  change  in  the  RMS 
variation  of  the  index  of  refraction,  O',  were  temperature  and 
salinity  microstructure,  and  bubbles. 
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The  effect  of  resonant  bubbles  on  the  change  in  a  would 
be  very  great.  Since  bubble  populations  decrease  with  depth,  this 
might  account  for  the  implied  decrease  in  01  with  depth.  When  in 
the  presence  of  a  motor  launch,  the  size  of  the  inhomogeneous 
volumes  were  made  smaller,  there  was  a  concurrent  increase  in  the 
bubble  population  and  in  the  value  of  Ct.  The  changes  and  vari¬ 
ations  in  the  other  parameters,  chiefly  temperature  and  salinity, 
should  be  accounted  for  in  subsequent  analyses,  and  a  quantitative 
estimate  of  the  relative  effect  of  bubbles  on  amplitude  modulation 
may  be  possible. 

5.  The  power  spectral  densities  indicate  that  most  of  the 
modulation  takes  place  at  those  frequencies  corresponding  to  the 
frequencies  of  the  surface  waves.  This  effect  was  noticable  down 
to  the  greatest  depth  (14.3  meters)  althouqh  the  prominemce  of 
the  PSD's  at  the  surface-wave  frequencies  tended  to  disappear  with 
greater  depth. 

6.  in  summary,  the  major  causes  of  acoustic  amplitude  modu¬ 
lation  at  60  kHz  in  this  near-surface  volume  of  ocean  could  well 
have  been  the  effects  of  mobile,  random  bubble  populations  acted 
upon  by  particle  velocities  generated  by  the  action  of  surface 
waves . 

F.  SUGGESTIONS  FOR  FURTHER  EXPERIMENTS 

The  basic  questions  left  incompletely  answered  by  this 
investigation  are  the  relative  effect  of  bubbles  versus  the  other 
causes  of  amplitude  modulation  near  the  surface,  and  the  spatial 
effect  of  the  inhomogcn iet ies  upon  the  transmission  of  sound. 

Later  analysis  of  the  oceanographic  data  will  help  to  answer  this 
question . 
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It  is  suggested  that  further  experiments  of  this  kind  would  be 
valuable.  In  order  to  investigate  the  spatial  effects  of  the 
inhomogen iet ies ,  the  basic  apparatus  of  this  experiment  could  be 
modified  as  shown  in  Figure  V-l6  in  order  to  attempt  some  cross- 
correlation  of  the  amplitude  modulation  received  by  each  of  the 
hydr  phones .  In  order  to  more  completely  investigate  the  question 
of  transport  versus  the  stationary  oscillation  of  inhomogenieties , 
more  thermistors  could  be  mounted  as  shown  on  Figure  V-l6. 

The  question  of  bubbles  and  bubble  populations  in  an  area 
might  be  answered  either  by  analysis  of  all  other  parameters  and 
their  effect,  or  by  deliberately  generating  bubbles  in  the  volume 
of  water  to  note  their  effect  on  ambient  conditions. 
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FOR  SPATIAL  STUDY 
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APPENDIX 

THE  RELATION  BETWEEN  F  AND  CV 


The  relation  between  F,  the  fractional  modulation  which  is 
used  in  the  development  of  the  experiment,  and  CV,  which  Stone 
and  Mintzer  call  the  "coefficient  of  variation"  will  be  developed. 

Fractional  modulation  is  defined  in  terms  of  the  voltages  as 
shown  below: 

p  -  fvMAX  -  V  ] 

F  l  VMAX  J 


where  V  =  modulated  voltage 
VMAX  =  unmodulated  voltage 

The  coefficient  of  modulation  is  defined  by  Stone  and  Mintzer 
in  terms  of  pressure  amplitude,  P,  as  shown  below: 

tcv)z  -  {<£>  -<p?\) 

<p>* 


Then  if  one  assumes  that  P  =  KV,  or  that  a  constant  factor 
relates  the  two  values,  as  it  does  in  the  equipment  used  for  these 
experiments,  then: 


(Cv) 


<vl> 


i_  _  £  [  <V;>  1 

<KV>1  **  I-  . 


Fractional  modulation  was  used  in  this  experiment  because 

small  fluctuations  were  anticipated.  A  development  between  the 

2 

variance  of  the  fractional  modulation  and  the  term  (CV)  is  shown 
below : 


F  =  Wuftx-vl  _  V_  \ 
L  VuiA*  J  '  ^Mflx  ) 
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[<F'>  '  <F>"] 

*Lf  (L.f?  -(<>■>-  <$..'>)" 

<<&.)*>  -  <A,f 


or,  since  V  is  a  constant, 

I’lnA 


6% 


(V  ikin') 


r<v*?  -  <v>* 


] 


<v> 


multiplied  by  the  factor 
expressed: 

<v>‘ 


<v>’- 


the  variance  may  now  be 


<V*>  -  <v>% 

""  <v>r~  J 


Where  the  second  term  is  a  numerical  value  multiplied  by 
Mintzer  and  Stones'  expression  for  the  coefficient  of  variation. 
The  value  multiplied  against  the  coefficient  of  variation  is 

<v>* 

which  is  numerically  obtainable  as  follows: 

<"?  >  -  1  -  <v> 


<V>  =  1  -  <F> 

for  all  runs,  the  value  of  V  was  6  volts. 

MAX 

0>l  =  (1  -  <f»x’ 

3C# 
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2  2 

The  relation  between  (T  and  (CV)  will  be  calculated  for  the 
4.3  meter  run: 

_  (I-  0.0SA6^_ 


-  0.0WH 


and  in  the  7.3  meter  run  was: 

<V>*  (1  —  _  O.C>2t£“3C> 

It  is  noted  that  the  conversion  factor  does  not  significantly 

change  with  depth,  and  therefore  conclusions  about  the  relative 
2 

value  of  (CV)  at  two  different  depths  are  applicable  as  well  as 

2 

the  relative  value  of  CT  at  these  depths. 


74 


nnnnnnno  ooooo 


PROGRAM  1 

PROGRAM  CONVERT!  INPUT, OUTPUT  .PUNCH) 

A  VGX  IS  THE  AVERAGE  DISTANCE  DISTANCE  THE  RECORDING 
DRUM  ADVANCES  PER  ONE  HOUR  REAL  TIME,  I.F.  THE  DATA 
SAMPLING  RATE, 

DIMENSION  UI6000) , V  (6000 ) , N ( 80  ) , I  BUFF ( 5000 ) , NK ( 80 ) 

IK  =  0 

READ  97  *  L 
97  FORMAT (12) 

READ  96.AVGX 
90  FORMAT ( F 1 0 . 0 ) 

READ  96,IFILE 
96  FORMAT (12) 

DO  116  J J= 1 , L 
DELT=36.0/AVGX 
PRINT  99 , OtLT 

99  FORMAT ( 1H0, 2X,25H  SAMPLING  INTERVAL  EQUALS ,  1X,E15.7, 
11X.7HSEC0NDS) 

PRINT  700, JJ 

700  FORMAT! 1H0.9HJJ  EQUALS,  2X,I2) 

DO  100  1=1,5500 
U( I ) =0 . 0 

100  V! I  1=0.0 

101  DO  202  1=1,5000 
202  I  BUFF  C I ) =0 

CUUNTX=0 .0 
COUNT Y=0 .0 
NUM=  5000 
M=  1 
K8  =  0 

CALL  L I  OF ( 5LRB  CD  1 , I BUF F , NUM , NP AR , NEOF ) 

IF(NEOF)  602,200 

200  K=— 7 
KF=0 

201  K=K+8 
KA=K+8 
KC  =  0 

DO  102  KB=K , KA 

kf-kf+l 

NK(KC)=IBUFF(KB) 

102  IF  ( NK(KCI .EO.O)  KF=1 

DECODE (80, 103, NK)  ( N( I ) , 1= 1 , 80 ) 

103  FORMAT ( 80R 1 ) 

DO  1 04  19=1,80 

IF(N< 19)  .EQ.50B)  GO  TO  106 
I F ( N ( 19) .EQ.55B)  GO  TO  107 
I F ( N ( 19) .EQ.34B)  GO  TO  108 

SYMCOL  /  (50)  REPRESENTS  AN  INCREMENT  TRAVEL  IN  THE 
MINUS  X  OR  Y  DIRECTION  BY  THE  DIGITIZER. 

SYMBOL  0, ( 55B ) ,  REPRESENTS  A  ZERO  INCREMENT  TRAVEL 

IN  THE  X  OR  Y  DIRECTION  BY  THE  DIGITIZER 

SYMBOL  1 , ( 34B ) ,  REPRESENTS  AND  INCREMENT  TRAVEL  IN 

THE  POSITIVE  X  OR  Y  DIRECTION  BY  THE  DIGITIZER 

SYMBOL  *,(47B),  IS  A  FLAG  INSERTED  IN  IN  THE  RECORD 

BY  THE  PERSON  DIGITIZING  BY  USE  OF  THE  IRG 

GO  TO  104 

106  RX=—0 .01 
K8=K  8  +  1 
GO  TO  109 

107  RX=0 . 0 
K8=K  8+ 1 
GO  TO  109 

108  RX=0 .01 
K8=K  8+ 1 

109  K3=K 6/ 2 
K3  =  2*K3 

IF(K3.EQ.K3)  GO  TO  111 
COUNTX=COUNTX*RX 
IF(COUNTX.NE.O.O)  GO  TO  110 
GO  TO  104 

110  U ( M ) =CCUNTX 
GO  TO  104 


111  COUNTY  =  CCJUNTY+RX 

IF (COUNT X. NE .0 . 0  )  GO  TO  112 
GO  TO  104 

112  V(M)=COUNTY 
C0UNTX=0 .0 
M=M+1 

IF(M.GT.bOOO)  GO  TO  600 

104  CONTINUE 
IF(KF.EQ.l)  GO  TO  113 
GO  TO  201 

113  MAX=M-7 

105  PRINT  205,  M,  19 
205  FORMAT! 1H0, 2110) 

C  TOTAL  TIME  OF  THE  RECORD  EQUALS  NUMBER  OF  DATA 
C  POINTS  TIMES  THE  SAMPLING  INTERVAL,  DELT. 
TIME=(M*DELT)/3600.0 
PRINT  115, TIME 

115  FORMAT ( 1H, 10X , 26H  TOTAL  TIME  OF  RECORD  EQUALS, 2X, 
1E15.7, 2X,7H  HOURS. ) 

PRINT  215, <V< I) , 1=1, M) 

PRINT  215, (U( I ), 1=1, M) 

215  FORMAT ( 1H  ,10X,14F7.2) 

PUNCH  213, ( Vt I ) , 1=1, M) 

213  FORMAT (14F5. 2) 

PUNCH  2100 

2100  FORMAT  ( 20H*#*******#*#****#**# ) 

116  CONTINUE 
STOP 

600  PRINT  601 

601  FORMAT ( 1H0.20X ,37H*****  U  AND  V  SPACE  INADEQUATE  ***,) 

602  IK= IK+ 1 

I F ( IK.GT.IFILE )603,200 

603  STOP 
END 
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PROGRAM  IV 

PROGRAM  TO  CALCULATE  THE  VAR 
MCDULAT I CN  OF  EACH  RUN. 
DIMENSION  X ( 56 00  ) 

READ  ( 5, 1100) ( X(  I  >  ,  1=1 ,5600) 
1100  FORMAT ( 14F5. 2 ) 

XMAX=X ( 1 ) 

DO  23  1=2,5600 
IF  IXMAX.GE.Xt I ) )  GO  TO  23 
XMAX=X ( I ) 

23  CONTINUE 

DO  24  1=1,5600 

X(n  =  (Xm-XMAX)*<-.212) 

24  CONTINUE 
X3  =  0 . 0 

DO  30  1=1,5600 
X3=X3+( X( I )**2  ) 

30  CONTINUE 
X3=X3/5600 

WR I T  E ( 6 , 99 ) ( X3  ) 

99  FORMAT ( 5X, *  ME  AN  SQUARE=',F5. 
X4=0 . 0 

DO  31  1=1,5600 
X4=X4+X { I ) 

31  CONTINUE 
X4=X4/5600 
WRITE ( 6 , 100) ( X  4 ) 

100  FORMAT (5X,  •MEAN=* , F5. 4) 
X5=X4**2 

X6=X3-X5 

WRITE(6,101)(X6) 

101  FORMAT ( 5X,  ' S IG MA  SQUARE=',F5 
STOP 

END 


VARIANCE 


FRACT ITNAL 


SQUARE= ' , F  5 .4 ) 


SQUARE= 1 , F5.4) 
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